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The distribution of the gluon action density in mesonic systems is investigated at finite temper- 
ature. The simulations are performed in quenched QCD for two temperatures below the deconfin- 
ment phase. Unlike the gluonic profiles displayed at T = 0, the action density iso-surfaces display a 
prolate-spheroid like shape. The curved width profile of the flux-tube is found to be consistent with 
the prediction of the free Bosonic string model at large distances. 



o ■ 
(N : 



i3 ■ 

Oh' 



> 
oo 

O 

^' 

O 

o 



X 
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I. INTRODUCTION 

In the flux-tube model, the linearly rising potential 
between a pair of color sources is believed to be due to 
the formation of a thin gluonic flux-tube of a constant 
cross-section. At high temperatures, lattice simulations 
on pure SU(3) gauge fields ^1^ have indicted a decrease 
of the effective string tension for the quark anti-quark 
potential with the rise of the temperature. The QCD 
vacuum structure around the sources is then expected to 
exhibit gluonic profiles with widths variant to the zero 
temperature case. The detailed geometry of the gluonic 
field at finite-temperature and whether it holds the con- 
stant cross-section property is an interesting topic that 
has not yet been explored in Lattice Quantum Chromo- 
dynamics (LQCD). 

The low-energy dynamics of the flux-tubes in the in- 
frared region of a confining gauge theory can be described 
in terms of an effective Bosonic string. The thin flux-tube 
between two widely-separated color sources fluctuates 
like a massless string. The linearly rising part of the po- 
tential arises from the classical configuration which corre- 
sponds to the flat world sheet of the string. The quantum 
fluctuations of the string lead to an universal sub-leading 
correction to the potential well known as the Liischer 
term 0]. Lattice simulations for several gauge theories 
[^-0] have supported the existence and the universality of 
the string's sub-leading effect. At high temperatures, the 
gluonic modes come into play and the effective string de- 
scription of the temperature dependent quark anti-quark 
potential has been worked out in Refs. [1, 0. A fur- 
ther comparison with SU(3) Monte Carlo lattice data for 
temperatures beginning from T = 0.8 Tc [l| has shown 
a good parameterizing behaviour to the string picture 
formula for a minimal distance of RT = 0.5. On the 
other hand, there has been numerical indications that 
the inclusion of the higher-order string effects beyond the 
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Gaussian approximation, e. g, the string's self-interaction 
terms in Nambu-Goto effective string action, has repro- 
duced the correct temperature-dependent string tensions 
up to a temperature scale of T = 0.5 Tc in the 3D gauge 
Z2 model [13. 

The string model predictions of a logarithmic broaden- 
ing [m for the width of the string delocalization has also 
been observed in several lattice simulations correspond- 
ing to different gauge groups 0, SHI, 111- As the temper- 
ature increases, substantial deviations from the logarith- 
mic behavior is expected, and the law of broadening turns 
eventually into a linear growth for large distances before 
the deconfinement is reached from below . Apart from 
the peculiar features in the laws of the broadening of the 
flux-tube when the temperature is raised, the string pic- 
ture predicts an effect which is rather interesting from 
the geometrical point of view. The width calculated at 
each corresponding transverse plane to the line joining 
the two quarks is found to differ from that at the cen- 
tral plane by an amount that increases with the rise of 
the temperature. In other words, the mesonic string pic- 
ture is implying a curvature in the gluonic proflle that 
becomes more pronounced as higher temperatures are 
approached. The string self-interaction with the quark 
line causes a noticeable difference in the delocalizations 
beyond the central transverse plane and these aspects re- 
main to be ascertained in LQCD. Moreover, revealing the 
whole profile of the glue at finite temperature provides 
a particularly interesting source of knowledge regarding 
the true geometry of the flux-tube, since at finite tem- 
perature one, naturally, need not hold to any particular 
assumption for the shape of the gluonic source wave func- 
tions in the relevant gauge-invariant objects representing 
the quark states. Probing the transverse profile of the 
glue might even be of relevance to the modeling of ground 
state sector of the theory where the exact geometry of the 
flux-tube seems to be not yet settled fT5|. 

In this paper, we investigate the distributions of the 
color fleld inside the meson at two temperatures below 
the deconfinement phase, T ~ 0.9 Tc, and T ~ 0.8 Tc. 
The lattice simulations are performed on the SU(3) 
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gauge group in the quenched approximation. The field 
strength inside the corresponding quark system will be 
revealed by correlating an improved action density op- 
erator [I6'] to the mesonic state. The mesonic states 
are accounted for by means of Polyakov-loop correlators. 
Gauge-smoothing [TtI, [I!] , in addition to a high statis- 
tics gauge-independent approach [l9| will be employed 
to enhance the signal to the noise in the flux correla- 
tion function. This noise reduction approach is variant 
to other approaches that utilize gauge- fixing [2^. The 
obtained profile of the action density will then be com- 
pared to the prediction of the mesonic string models at 
several distances for the highest temperature near the 
deconfinment point T ~ 0.9 Tc. 

The paper is organized as follows: In Sec. II, the details 
of the simulation will be described. We review the pre- 
dictions of the Bosonic string model for the qq potential 
and the width of the string fluctuations at finite tem- 
perature in Sec. III. In Sec. IV, we measure the quark 
anti-quark potential and examine the string model pa- 
rameterization for various levels of gauge smoothing. In 
Sec. V, the action density in mesons will be presented, 
the corresponding widths at several transverse planes to 
the tube is then measured and set in comparison with 
the string model predictions. Sec. VI provides the con- 
clusion. 



by a gauge-invariant action density operator S'(p, t) at 
spatial coordinate p of the three dimensional torus corre- 
sponding to an Euclidean time t. The measurements are 
repeated for each time slice and then averaged. 
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(2) 



nt = l 



The action density operator is calculated via a highly- 
improved 0{a*) three-loop improved lattice field strength 
tensor [l6j. A dimensionless scalar field that character- 
izes the gluonic field can be defined as. 



C(p;ri,r2) 



;^2Q(ri,r-2)) {S{p)) - {V2Q{ri,r2)S{p)] 
{V2Q{ri,r2)) {S{p)) 



(3) 

where < > denotes averaging over configurations 

and lattice symmetries, and the vector p refers to the 
spatial position of the fiux probe with respect to some 
origin. Cluster decomposition of the operators leads to 
C ^ away from the quarks. 



Noise reduction 



II. SIMULATION DETAILS 

A. Color field measurements 

In this investigation we have taken our measurements 
on 500 quenched QCD gauge-field configurations for each 
set of lattice parameters considered. The gauge configu- 
rations were generated using the standard Wilson gauge 
action on lattices with spatial volume of 36'^. We chose to 
perform our analysis with lattices as fine as a = 0.1 fm by 
adopting a coupling of value /3 = 6.00, with temporal ex- 
tents of Nt — 8, and Nt = 10 slices, which correspond to 
temperatures T ~ 0.9 Tc, and T ~ 0.8 Tc, respectively. 
The gluonic gauge configurations were generated with a 
pseudo-heatbath algorithm [U upda ting the corre- 
sponding three SU(2) subgroup elements [23. Each up- 
date step consists of one heat bath and 4 micro-canonical 
refiections. The measurements are taken after each 2000 
of updating sweeps. 

The mesonic state is constructed by means of Polyakov 
loop correlators, 

V2Q = {P{n)P\r2)), 
where the Polyakov loop is given by. 



Pin) = ^Tr 



Nt 



(1) 



and the vectors define the positions of the quarks. The 
measurements that characterize the color field are taken 



We make use of translational invariance by comput- 
ing the correlation on every node of the lattice, averag- 
ing the results over the volume of the three-dimensional 
torus, in addition to the averaging the action measure- 
ments taken at each time slice in Eq. ([2]). To further 
improve the signal to noise ratio in the gluonic corre- 
lation function, local action reduction by smearing the 
gauge links has been performed on the whole 4 dimen- 
sional lattice. Since the main focus in this investigation 
is to resolve the nature of the flux distributions in the 
IR region of the theory, we have been able to show that 
with an appropriate level of gauge smoothing, effects on 
the large distance correlations can be kept minor. For 
each distance scale, a range of smoothing levels is seen to 
leave physical observables and topological structures in- 
tact. Similar techniques have been adopted in Ref. in 
the determination of the large distance QQ force in vac- 
uum with different levels of HYP smearing. In Sec. IV- 
A, it has been shown that for a given distance scale, the 
measured quark anti-quark force at large distances can 
be left with negli gible changes for a range of smoothing 
levels. Variant to [2^ where the Cabbibo-Marinari cool- 
ing has been employed, we have chosen to smooth the 
gauge field by an over-improved stout-link smearing al- 
gorithm [1^. The use of this algorithm should ensure 
that the smoothing algorithm has a minimal effect on the 
topology of the gauge field [11]. In standard stout-link 
smearing [l3], all the links are simultaneously updated. 
Each sweep of update consists of a replacement of all the 
links by the smeared links 



Ufi{x) = exp(i(3^,(x)) Uf_t{x) , 



(4) 
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with 

-^tr(17t(x)-f]^(x)), 

and 

where Yj^^{x) denotes the sum of the two staples touching 
U^{x) which reside in the n — v plane. The scheme of 
over- improvement requires E^i,(a;) to be replaced by a 
combination of plaquette and rectangular staples. This 
ratio is tuned by the parameter e [18|. In the following 
we use a value of e = —0.25, with — p = 0.06. We 
note that for a value of p = 0.06 in the over-improved 
stout-link algorithm is roughly equivalent in terms of UV 
filtering, to the standard stout-link smearing algorithm 
with the same p = 0.06. 

III. THE BOSONIC STRING MODEL 

A. Temperature-dependent quark— ant i- quark 
potential 

The correlation function of two Polyakov loops on the 
lattice determines the interaction potential between the 
color sources. 



where gap is the two dimensional induced metric on the 
world sheet embedded in the background , 

9 = det(ga/3). 

Gauge fixing is required for the path integral (9) to be 
well defined with respect to Weyl and reparamterization 
invariance. The physical gauge, X^ = Ci,-^** = C2 would 
restrict the string fluctuations to transverse directions to 
C. In the Quantum level, Weyl invariance is broken in 4 
dimensions, however, the anomaly is known to vanish at 
large distances [2^. The action after gauge fixing read, 

S\X\=a r dC, rdC2(l + (%^i)' + (%^±)')i 
Jo Jo 

(8) 

Expanding the square root in powers of aRL, 

S[X]^aRL+^ TdCi r dQ2{yxf + , (9) 

the action decomposes into the classical configuration 
and fluctuation part, and the string higher-order self in- 
teractions. A leading order approximation can be made 
by neglecting the self-interaction terms, the path integral 
Eq. ^ is then. 



(P(O)pt(i?)) = j d[U]P{Q)P\R)eM-S^), 

= eM-V{R,T)/T). (5) 

Su] is the plaquette action and T is the physical tempera- 
ture. The self interactions of the glue exchanged between 
two color sources in QCD can result in the squeezing 
of the glue into a thin one dimensional string-like ob- 
ject. The immediate consequence of this string picture is 
that a functional form can be ascribed to the Polyakov- 
loop correlators. Namely, the partition function of the 
string. The correlators are expressed as a functional in- 
tegrals over all the world sheet configurations swept by 
the string, 

(P(O)pt(i?))^ / [DX]cM~S{X)). (6) 
Jc 

The vector X''(Ci,C2) maps the region C C R^ into 
R^, with Dirichlet boundary condition X{(iX2 = 0) = 
^(CiiC2 — R) = and periodic boundary condition 
along time direction X{(i = 0, C2 = 0) = X{(i = i, C2 = 
-R), and S is the string action and can be chosen to be 
proportional to the surface area, i.e. the Nambu-Goto 
action, 

S[X] [dCi [ dC2 V9, (7) 



(P(0) P\R)) = e-"«^ [det(--V2)]-i. (10) 

The determinant of the Laplacian on the cylinder has 
been regulated using a lattice regulator in [8]. The po- 
tential is obtained in closed form for a length scale com- 
parable to the thermodynamic scale in [9]. The effective 
potential is, 

V{R, T) = (^(j -^T^ +'^T^ tan"\2 i?T)"^^ R 

_ _itan-H2 i?T)^ - 
V12 6 ^ 7 i? 

- I ln{l + {2 RTf)+f,. (11) 

The limit of large string length entails taking the 
temperature-dependent string tension to be, 

a{T)^a-^T\ (12) 

The free string model predicts a temperature-dependent 
quark anti-quark potential that is featured by the exis- 
tence of a logarithmic term in addition to a leading-order 
decrease in the string tension by an amount ^ . 
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B. Width of the Flux-tube fluctuations at finite 
temperature 



with the nome, q = e 2 The expectation value of the 
mean square width would then read, 



The vibration modes of the string-like object renders an 
effective width for the flux-tube. A well known predic- 
tion made by Lusher, Miinster and Weisz based on 
the effective Bosonic string model, has shown that the 
mean square width of the vibrating flux-tube at the cen- 
ter plane grows logarithmically as a function of the in- 
terquark separation, R, 



w 



— log(T) 



(13) 



where A is an ultra-violet scale. With the increase of 
the temperature, higher order gluonic degrees of freedom 
are present and the effective width of the correspond- 
ing string is expected to manifest an intricate behavior 
involving both the distance and the temperature. The 
mean square width of the string is defined as, 



Jc[DX]X^eM-S[X]) 
L[DX]eM-S[X]) ' 



(14) 



C = (^11*^2) is a complex parametrization of the world 
sheet, such that G [-R/2,R/2],^2 & [-L/2, L/2], with 
T = is the modular parameter of the cylinder, and L is 
the temporal extent governing the inverse temperature. 
Casselle et al. [0] and Gliozi _27| have worked out the delo- 
calization of the string for all the planes transverse to the 
line joining the quark pair by the corresponding Green 
function. This technique proceeds by removing the diver- 
gence in the quadratic operator in Eq. p^ by the use of 
the Schwinger [28!| point-split regularization, then taking 
the limiting action for the Numbu-Goto model as that 
of the corresponding Gaussian model. The quadratic op- 
erator is then the correlator of the free Bosonic string 
theory in two dimensions. 



{X^^;r)) = {X{0-X{^ + e)), 

= G(e,e + e). 



(15) 



This Green function is the solution of Laplace equation 
on cylinder with a Dirichlet boundary condition. 



G(C,eo)=^log|/(C,eo)|. 

The conformal map reads 

9i[7r(e-6)/i?;T] 



/(e,eo) 



92[7r(e-eo)/^;^]' 
where the Jacobi 9 functions are, 



(16) 



(17) 



01 (C; r) = 2 ^(-l)"q"("+i)+i sin((2n + 1) 0, 

00 

92{(.; t) 2 ^q"("+i)+3 cos((2n + 1)0, (18) 



1 1 



27rCT 



log 



g2(7rgi/i?;T)) 
9[{0;t) 



(19) 

This expression converges for modular parameters close 
to 1, and contains in addition to the logarithmic diver- 
gence term a correction term that encodes the depen- 
dence of the width at different transverse planes on the 
modular parameter of the cylinder. At finite tempera- 
ture, this term is contributing to the width at all the 
planes. Fig. 1 is a plot of the mean square width calcu- 
lated at ^1 values via Eq. (fT9| . The plot shows the profile 
for several modular parameters and fixed separation be- 
tween the two Polyakov loops. 




FIG. 1: The mean square width, Ea. (fT9|) . of the flux-tube 
evaluated at all planes ^1 perpendicular to the quark-anti- 
quark line. The separation distance between the pair is i? 
= 10. 



The string model predicts an increase in the width with 
the increase of the temperature. The increase in the 
width is maximum at the central plane which is seen as an 
increase in the curvature in the profile of the string fluctu- 
ations. At zero temperature L — )■ 00, Eq. ([T^ converges 
well, and the second term in Eq. still contributes to 
the whole shape of the fluctuations at all planes except 
the middle, the contribution of this term at zero temper- 
ature is, 



1 1 I 

^log|cos( — ) 



(20) 



n=0 



which is seen from the plot Fig. 1 to be subtle in the 
middle region and have more pronounced effects on the 
width near the quark positions. 
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IV. RESULTS 
A. Quark— anti-quark potential 

The technique adopted to enhance the signal to the 
noise in the correlation function which characterizes the 
gluon flux, Eq.Q, involves smoothing the gauge links 
by the over-improved stout-link smearing algorithm de- 
scribed above. The whole 4-dim torus has been smeared 
for the consecutive levels of smearing correspond to 20, 
40, 60, and 80 sweeps, forming four data sets for in each 
temperature. The choice of the appropriate data set 
for the numerical evaluation of the expectation values in 
Eq. ^ at each distance scale, should be based on a com- 
promise to simultaneously achieve two tasks, the smear- 
ing level has a minimal effect on the physical observables, 
and a significant error reduction is gained. The larger the 
separation distance between the quark pair, the higher 
smearing level required to gain good signal to noise in the 
correlations in Eq.Q. However, smearing has an effect 
on the observables equivalent to the increase of the lattice 
space-time cut off, and a large enough number of smear- 
ing sweeps will result in a subsequent loss of the physics 
on the short distance scale. The physical observable of 
direct relevance to the properties of the gluonic flux-tube 
is the quark anti-quark potential. For each level of smear- 
ing, we numerically evaluate the quark anti-quark poten- 
tial and the corresponding force. At fixed temperature 
T, the Monte Carlo evaluation of quark-anti quark po- 
tential at each R is calculated through the Polyakov loop 
correlators according to. 



V{R,T) 



-1 



log((P(0)pt(i?))). 



(21) 



The jacknife error analysis for the data shows a sig- 
nificant decrease in the uncertainties associated with 
Polyakov-loop correlators on short distance scale when 
measurements are taken after 20 sweeps of smearing. 
For large distances, a subsequent increase of 20 sweeps 
would provide error reduction by factors of 1.3 < x < 1.5 
for the corresponding distances 0.6 fm < i? < 1 fm. 
Table I summarizes the factors of error reduction for the 
Polyakov loop correlator after each incremental increase 
of 20 sweeps of smearing. j39| 



TABLE I: The error reduction factor in the Polyakov-loop 
correlator Eq. ([5]) by the increase of the number of 20 smearing 
sweeps for each smeared data set. 



R /No. Sweeps 


20 - 40 


40-60 


60 - 80 


6 


1.3 


1.1 


1.1 


8 


1.4 


1.2 


1.2 


10 


1.5 


1.3 


1.2 



To test the validity of the gauge smoothing approach, 
or equivalently, to determine the levels of smoothing for 
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(b) 



FIG. 2: The quark-anti quark potential measured at each 
depicted smearing level, the lines correspond to fits of the 
potential obtained from the string picture of Eq. (|lll) for each 
data set as described in the text. The upper plot is at T = 
0.8 Tc while the lower plot is at T = 0.9 Tc. 



which the physics is left intact, one is tempted to set a ref- 
erence scale which signifies how the smeared data would 
behave with respect to the string model parametrization. 
This approach of referencing the data to the string model 
is justified by the fits previously reported in Ref.[]|, which 
has returned good and shown stability to the fit range 
at large distances. 

The numerical data obtained for the quark anti-quark po- 
tential Eq. (PT|) on every smoothed gauge configuration is 
fitted to the string picture qq potential Eq. (fTT|) . The 
effects of smearing is expected to be more pronounced at 
short distances. For this reason the minimal fit distances 
is taken as large as possible, R > 0.7 fm for T = 0.8 Tc, 
and R > 0.9 fm for T = 0.9 Tc. The string tension has 
been taken as a fitting parameter. The fits are return- 
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ing good for all the smoothed data sets considered. 
The quark anti-quark potential and the corresponding 
fits are shown in Fig. 2. The fits to the data show almost 
equal slopes for all smoothing levels. This is also mani- 
fest in Table II, where the string tensions are measured 
in accord to Eqs. (11) and (12). Within the standard 
deviations of the measurement, the string tensions for 
all levels of smearing are equal values. At temperature 
T = 0.9 Tc, our measurements for the string tension agree 
for all the data sets. Moreover, this value is in agreement 
with that reported in Ref . [l| . The factors of error reduc- 
tion at higher temperature at i? = 1 fm after 40 sweeps 
of smearing equals the corresponding one at T = 0.8 
after 80 sweeps. The noise tends to decrease with the 
increase of the temperature. This analysis shows that 
for the qq separation distances depicted in Table II, all 
the smoothed configurations are appropriate for revealing 
the gluonic field. On the other hand, the data points for 
R < 0.7 at T = 0.8 shift up towards the string model's 
curve with the increase of the number of smearing sweeps. 
The removal of short distance physics is manifest here. 
The difference in the regularization brought about by the 
increasing of the space time cut-off introduced by smear- 
ing, shifts the qq potential by a renormalization constant 
in Eq. (pij) . To manifest the effect of smearing on the 
qq potential, the potential Eq. ([21]) has not been nor- 
malized. The qq force, however, can be calculated to 
eliminate these constant shifts. With the definition of 
the derivative on the lattice taken as Refs. [3, H^, the 
force is computed as. 



V{r) ~ V{r - a) 



(22) 



Fig. 3 shows the force calculated for all smearing levels 
for distances up to 1.4 fm. The force from the string pic- 
ture, Eq. ((TT|) , with fit parameters measured at 60 sweeps 
in Table II. is illustrated. The qq force calculated from 
the string model compares well with that data measured 
on 20 sweeps of smearing for distances commencing at 
0.55 fm. All data points correspond to smearing levels of 
20, 40 and 60 sweeps coincide, within the statistical de- 
viations, for distances greater than 0.75 fm. For 20 and 
40 sweeps the data convergence begins one lattice spac- 
ing earlier at R=0.65 fm. The essential features of the 
confinement remain unchanged for the considered levels 
of gauge smoothing on the scale of large distances. 
A combination of a large number of Monte Carlo updates 
followed by link averaging has been performed in Ref. [l[ 
to evaluate the qq potential for a range of temperatures 
above and below the confinement phase. This involves 
a large number of updating sweeps and measurements 
which makes it rather expensive in terms of the CPU 
time. This is particularly true for the evaluation of the 
gluonic flux distribution, since this concerns the Monte 
Carlo evaluation of not only the Polyakov-loop correla- 
tor, but also the three-point correlation function in the 
numerator of Eq.(|31). Gauge smoothing is chosen as a 



TABLE II: The string tension measured on all data sets cor- 
responding to various levels of link smearing. The measure- 
ments are obtained from the fits to Eq. (11) and (12). 



No. Sweeps 




Fit Range n — R/a 


T = 0.8 Tc 






20 


0.047(3) 


8-12 


40 


0.050(2) 


8-12 


60 


0.0493(9) 


8-12 


80 


0.0478(6) 


8-12 


T = 0.9 Tc 






40 


0.0385(8) 


10-13 


60 


0.0377(9) 


10-13 


80 


0.0373(8) 


10-13 



0.18 
0.16 
0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 





Model 
SW20 
SW40 
SW60 
SW80 



i— I- 



0.2 



0.4 



0.6 O.J 

R(fm) 



FIG. 3: The qq force measured for all the smearing levels up 
to distance of 1.4 fm, the temperature is T = 0.8Tc , /3 = 6. 
The line denotes the force as predicted by the string model 
at finite temperature, Eq. (|11[) . 



cheap and effective method in this case to reveal the gen- 
eral topological features of the flux distribution which can 
be confronted with the predictions of the string model. 
We have been able to show in this section the ranges of 
the validity of this approach, through the measurements 
of the physical observables that have been previously re- 
ported in [H,i.e the qq potentials and the string tension. 



B. Action Density 

1. Tube profile (qualitative picture) 

The lattice operator which characterizes the gluonic field 
is usually taken as the correlation between the vacuum 
action-density S{p,t), and a gauge- invariant operator 
representing the quark states. At finite temperature this 
must be a pair of Polyakov lines. The action-density 



(a) R = 5 a 



(h) R = 6 a 





(c) R = 7 a 



(d) R = 8 



FIG. 4: The flux distribution C{p,r\,f2) as given by the characterization Eq. ((3)) in the plane of the quark anti-quark pair 
p{x,y,z — zq), for separation distances R (a) 0.5 fm, (b) 0.6 fm, (c) 0.7 fm and (d) 0.8 fm at T = 0.8 Tq. The spheres refer 
to the position of the quark and anti-quark. 



operator is calculated through an O(a^) improved lat- 
tice version of the continuum field-strength tensor. Dis- 
cretization errors are reduced by combining several clover 
terms complemented by tadpole improvement [l^. We 
take our measurements with a three-loop field-strength 
tensor given by, 



(23) 



where C'*^*' is a combination of Wilson loop terms cor- 
responding to loops with lattice extent i used to con- 
struct the clover term and Wi are weights [l6j . The re- 



constructed action density, 



(24) 



is accordingly measured on 20 sweeps of stout link smear- 
ing. This has the effect of the removal of the divergence in 
the action density in the neighborhood of the quark posi- 
tions. It is, however, very beneficial in obtaining a good 
signal to noise to display the flux strength. The correla- 
tion function Eq. ^ is found C(/5) > 0, and C ~ away 
from the quark position. The scaled flux distribution in 
the plane of the qq pair is plotted in Fig. 3, for several 
qq separation distances R = | ri — r2 |, at temperature 
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(a) R = 9 a 
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» 



(b) R = 10 a 

FIG. 5: The flux contour-line distribution in the plane of 
the quark anti-quark pair zq, for separation distances of 
(a) 0.9 fm, (b) 1.0 fm, the spheres denote the position of 
the qq pair, T = 0.9 Tc. 



T = 0.8 Tc. The Polyakov-loop correlator is measured 
on 40 sweeps of link smoothing. The distribution shows 
a peak in the middle point between the qq pair at small 
separation distance R = 0.5 fm. As the two quarks are 
pulled apart, the distribution C{p) decreases rapidly, the 
peak behavior diminishes and the distribution is almost 
constant at i? = 0.8 fm. The qualitative description of 
these density plots suggests a two-dimensional Gaussian- 
Uke behavior, however, as we will see in the next section, 
careful measurements of the widths at each perpendicu- 
lar plane to the qq line, yield different widths for large 
distances. 

The behavior of the flux-distribution around the outer 
edges of the density profile does depend on finite- volume 
[30|. As a by product of performing the simulations 
on large lattice sizes to gain high statistics in a gauge- 




FIG. 6: The flux iso-surface passing through the quarks, plot- 
ted together with a surface plot for the density distribution 
in the qq plane (Inverted). The measurements are taken on 
80 sweeps of smearing for separation distance R = 9a, and 
T = 0.8 Tc. The lattice spatial extent is 36^, /3 = 6. 



independent manner, the two lattices employed in this 
investigations are being of a typical spatial size of 3.6'^fm'^ 
which do minimize the volume effect. 
The curvature in the flux lines is manifesting itself as it is 
evident from the flux-contour plots in Fig. [51 The contour 
plots reveal the form of the flux tube just before the de- 
confinement phase T = 0.9 Tc, for qq sources separated 
by i? = 0.9 fm and R = 1 fm, respectively. A simi- 
lar plot of the action-density iso-surface at i? = 0.9 fm 
in Fig. [6] displays a three-dimensional version of Fig. [5] 
(prolate-spheroid like shape) for the flux-tube. This geo- 
metrical form of the density plot manifests itself at tem- 
perature T = 0.8 Tc which is known to be near the end 
of the plateau of the QCD-phase diagram [31i] . 
It is worth noting, nevertheless, that at zero temperature, 
the correlation of the action density with the Wilson loop 
taken as a mesonic operator, does not reveal this curva- 
ture of the flux lines in the inner region between the qq 
pair at large separation distance p|. Thus we have illus- 
trated how thermal effects show up in the action-density 
correlations for the first time to the best of our knowledge 

2. Tube profile (quantitative aspects) 

Usually studies carried out on the flux-tube laws of 
growth focus their measurements on the central plane 



9 



transverse to the qq line. At T = 0, it seems also that 
there is a wide belief that the tube has almost constant 
cross-section with a uniform energy- density profile for 
large qq separations. Nevertheless, at high temperature 
where the string tension is reported to decrease by a value 
around 10% at T = 0.8 Tc 1], our calculations of the 
flux chromo-strength inside the meson Figs. 4-6 display 
a non-uniform action-density pattern around the whole 
qq line. 

It has been conceived a long time ago [s^-fs^l that the 
QCD vacuum behaves like a dual superconductor, and 
the color field generated by a pair of quark sources is 
squeezed into a thin string-like object dual to Abrikosov 
vortex by the dual Meissner effect. This squeezed flux- 
tube has been conjectured [2] that it can vibrate as a free 
string. At high temperature, one would expect higher 
modes relevant to the collective degrees of freedom of the 
string-like object to give rise to new interesting measur- 
able effects, which seem not only to be related to the law 
of the growth of the tube's width IJ] , but also the width's 
proflle it self. The string model's solution, Eq. (|19l) . in- 
forms us about how the tube would behave behind its 
symmetry point in the middle, together with the observed 
chromo-fleld profile. Figs. 4-6, one is tempted to investi- 
gate this string effect and establish a quantitative com- 
parison between the model and the glue profile in QCD. 
This is the aim of this section. 

Different possible components of the field-strength tensor 
in Eq. (|23p can separately measure the chromo-electric 
and magnetic components of the flux. The action density, 
however, is related to the chromofields via ^{E'^ — B^) 
and is the quantity of direct relevance to the comparison 
with the string fluctuations Eq. ([13]). The width of the 
flux-tube may be then estimated through fitting the den- 
sity distribution C{p), Eq. ([3]), in each selected transverse 
plane p{xi,y,zo) to a Gaussian l6l, [T]|. 
The width of the tube is defined as, 



J (P rjrj 



(25) 



77 is the set of vectors perpendicular to the qq line in the 
Xi plane. 

The flux calculations with Polyakov lines as a mesonic 
operator are well known to be distorted by statistical 
noise. To take reliable measurements to reveal the tube's 
fine structure, we choose to perform our analysis on the 
tube's width at the highest temperature T — 0.9 Tc 
where the scalar field C{p) has smaller jackknife error 
bars, even at very large distances. To further suppress 
the statistical fluctuations, the number of measurements 
has been increased by a factor of 4. This has been done 
by updating the raw gauge flelds and then repeating the 
measurements described in Sec. II after each 70 sweeps 
of Monte-Carlo updates. The density distributions have 
been symmeterized around all the symmetry planes of 
the tube, the resultant average density C(p) is fit to a 
Gaussian of the form A{xi) e"'^"^") , with yo on the 
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FIG. 7: The density distribution C{p} for separation distance 
of _R — 12a, T = 0.9 Tc, plotted for the transverse planes 
X — 1, X — 3, and x — 6. The lines correspond to the Gaussian 
fits to the density in each plane p{xi,y, zq). 



qq line, see e.g Fig. 7. The Gaussian fits to the data are 
for several transverse planes between two sources sepa- 
rated by a distance of i? = 12 a. Table lllll summarizes the 
measurements on both the widths w'^{xi) and the ampli- 
tudes A{xi) of the fiux tube, in accord to these Gaussian 
fits at each transverse plane Xi to the qq line. The co- 
ordinates lattice co-ordinates (lattice units) and 
are measured from the quark position a; = 0. The uncer- 
tainties in width measurements at each transverse plane 
are the standard asymptotic errors in the Gaussian fits 
and can be assumed to be correlated. The fiux-density 
measurements at each source separation are all taken on 
40 sweeps of smearing. According to the discussion in 
Sec. IV-A, this level of gauge smoothing should leave 
the qq potential and force with insignificant effects for 
R > 0.6 fm. 

For a fixed source separation, the measured values in Ta- 
ble [ITT] and [IV] are indicating, generally speaking, changes 
in the tube width along the qq line. The maximum width 
is measured at the tube's symmetry point in the middle. 
At relatively small separations R < 0.9 fm, the change 
in tube width along the planes is subtle. It is seen to 
vanish at R = 0.6 fm and R = 0.7 fm. The variation in 
the tube's width, however, is more pronounced at large 
source separation distances, see e.g Table IIVI in general 
qualitative agreement with the predictions of the string 
picture Fig. 1. The growth in width with increasing the 
source separation is also maximum at the tube's center 
point. 

Since we focus here on comparing the tube geometry to 
the string profile rather than the laws of growth, we cir- 
cumvent any ambiguity in measuring model's fit param- 
eters Eq. (|19p . by measuring the change in the width of 
the tube at each corresponding plane with respect to the 
central plane xq, 
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TABLE III: The width and the amphtude of the flux-tube at each consecutive transverse planes Xi from the quark to the 
middle of the qq line. The measurements for sources separation distances R = 6 a to R = 10 a, for the temperature T — 0.9 T^. 



plane 
n=R/a 


X = 

A 


1 


X 

A 


= 2 

2 —2 

w a 


X 

A 


= 3 


X 

A 


= 4 


6 


0.093(1) 


15.6(4) 


0.108(1) 


15.6(4) 


0.113(1) 


15.7(3) 






7 


0.099(1) 


16.6(4) 


0.116(1) 


16.6(4) 


0.125(1) 


16.8(3) 






8 


0.101(2) 


17.3(6) 


0.120(2) 


17.6(5) 


0.131(2) 


17.9(5) 


0.135(2) 


18.0(5) 
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0.102(2) 


18.1(6) 


0.120(2) 


18.6(6) 


0.132(2) 


19.1(5) 


0.138(2) 


19.3(5) 
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(b) _R = 8 a 
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(c) _R = 9 a 



(d) R = 10 a 



FIG. 8: The width difference Sw'^ — w'^{xi) — w'^{xo) for qq separations (a) 0.7 fm, (b) 0.8 fm, (c) 0.9 fm, and (d) 1 fm, 
, 13 — 6, T = 0.9 Tc- The line denotes the width difference Sw^ as predicted by the string model Eq. (I19|l . 



6w^ =w^{x^)-w^{xo), (26) 

this can provide a measure on how rounded or squeezed 
the flux tube would be compared to the width of the 
string fluctuations. Fig. 8 shows the change in the 
tube width calculated for separation distance R = 0.7 
fm to i? = 0.9, with uncertainties taken assuming the 
standard errors in the Gaussian fits are correlated, i.e. 
I e{xi) — e{xQ) \. In contrast with the predictions of 
the string model, the tube has almost constant width at 
R = 0.7 a, the measured changes in width at the plane 



X = 1 deviate from the model at i? = 0.8 fm and R = 0.9 
fm by large values of 82% and 70%, respectively. The 
deviations decrease as the sources are pulled apart to 
54% at R = 10, 38% at i? = 12, and good agreement 
between both profiles is reached at i? = 14 as can be 
seen in Fig. 9. The change of the width measured at the 
inner-transverse planes, however, agrees with the model 
at shorter distances, i? = 12 for the plane x — 2 and 
a; = 3. In general, the four plots in Fig. 9 show significant 
improvement with respect to the model predictions com- 
pared to the four plots at shorter distances in Fig. 8. The 
fiux-tube shows a constant cross section for R = 0.6, 0.7 
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(b) R = 12 a 
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FIG. 9: Similar to Fig. 8, the change in width is plotted for qq separations (a) 1.1 fm, (b) 1.2 fm, (c) 1.3 fm, and (d) 1.4 fm. 



TABLE IV: Similar to Table. IIIIl the widths of the flux-tube are measured at each consecutive transverse planes Xi from the 
quark to the middle of the qq line. The measurements for sources separation distances i? = 10 a to 7? = 13 a. 



plane 


X = 1 


a; = 2 


X = 3 


X — A 


a; = 5 


a; = 6 


n=R/a 


2 -2 

w a 


2 —2 

w a 


2 —2 

w a 


2 -2 

w a 


2 —2 

w a 


2 -2 

w a 


10 


18.7(7) 


19.4(6) 


20.1(6) 


20.5(6) 


20.6(6) 




11 


19.0(5) 


20.0(6) 


21.0(6) 


21.6(6) 


21.9(6) 




12 


19.2(7) 


20.3(6) 


21.6(6) 


22.6(7) 


23.2(6) 


23.3(6) 


13 


19.3(6) 


20.3(5) 


21.7(7) 


23.1(7) 


24.1(7) 


24.6(7) 



fm in disagreement with he string picture. At distances 
0.7 f ni < R < 1.0 fm, the LGT gluonic distribution profile 
is, geometrically speaking, more squeezed than the free- 
string picture would imply. As the sources are pulled 
further apart, the disagreement decreases gradually and 
the profiles of the glue and the string both compare well 
for sources separations i? > 13 fm. 

The thermal effects are manifest in the gluonic profile, 
giving rise to non- uniform widths. The string picture can 
parameterize these profiles only at large distances. At 
short distances on the other hand, the free-string picture 
does not seem to model the gluonic interactions on the 
scale of short distances which may become even more 



relevant in the thermal regime. 



3. Tube growth in width 

The measured values in Tables IIIIl and IIVI indicate a 
growth in the tube's mean square width at all trans- 
verse planes Xi as the color sources are pulled apart. 
The growth in flux-tube width at each selected trans- 
verse plane can be compared to the corresponding growth 
in the string fluctuation Eq. (|19p . this comparison can 
be performed by the fitting of formula Eq. ([T5| to the 
tube measured widths. Table IVl summarizes the resultant 
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{h)x: 
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(C) X : 
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FIG. 10: The width w'^{xi) for qq separations i? = 6 a to i? = 13 a at four consecutive planes (a) x = 1, (b) x — 2, (c) x = 3, 
and (d) x = 4. P — 6 , T — 0.9 Tc . The line denotes the string model, Eq. (|19p . fit to the data as described in the text. 



TABLE V: The resultant measurements of the scale Ro for the first four consecutive transverse planes Xi in accord to the fits 
of the tube width to the string model formula Eq. 1191 The values of the fit parameter and the corresponding Xdof a-re presented 
for variety of fit ranges. 



Plane 


x^l 




X — 


2 


a:: = 3 




X ^ 4 




Fit range 


Ro 


2 

Xdof 


Ro 


Xdof 


Ro 


Xdof 


Ro 


Xdof 


6-10 


0.0155(3) 


4.3 


0.028(5) 


5.9 


0.033(5) 


5.5 






6-13 


0.0126(2) 


5.9 


0.022(4) 


8.9 


0.027(4) 


9.2 






7-13 


0.0108(1) 


3.1 


0.019(3) 


5.7 


0.023(4) 


7.1 






8-13 


0.0093(8) 


1.2 


0.017(2) 


2.4 


0.019(3) 


3.1 


0.020(3) 


4.7 


9-13 


0.0086(5) 


0.4 


0.014(1) 


0.9 


0.017(2) 


1.7 


0.018(3) 


2.9 


10-13 


0.0081(2) 


0.1 


0.0129(6) 


0.3 


0.015(1) 


0.7 


0.015(1) 


1.2 



measurements of the fit parameter and the correspond- 
ing Xdof ^'^^^ consecutive transverse planes a; = 1 to 
x = 4. The fits show strong dependency on the fit range if 
the points at small sources separations are included. The 
highest value of is returned when fits include the whole 
range of sources separations, i.e. i?=:6toi? = 13. With 
the first four points excluded from the fit, the returned 



Xdof smaller indicating that the data points at large 
source separation are better parameterized by the string 
model formula. The value of the Xdoi gradually decreases 
as we exclude points at short distance separations, and 
stability in the fit is reached for widths measured for the 
plane a; = 1 at sources separations R > 0.7 fm, and at 
R > 0.8 fm for the plane x = 2. The fits are returning 
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good Xdof values for fits at the planes a; = 3, 4 for sources 
separations R > 0.9 fm. In the regions where the fits 
are returning good Xdof ' ^he values of the fit parameters 
are almost equal for the planes in the middle, a: = 2, 3, 4. 
However, at the closet plane to the sources, x — 1, the 
value of the returned parameter, unsurprisingly, deviates 
from the corresponding one at other planes. This is a 
manifestation for the above observed deviations in the 
change in tube widths at this plane compared to the cen- 
tral plane Figs. 8 and 9. 

Fig. 10 shows data points and the corresponding best 
fits to the string model at each plane, the string model 
at finite temperature poorly describes the lattice data at 
short distances. The plots depict the fact that the flux 
tube observed in LGT has a more suppressed profile than 
the fiuctuations of the free string would imply at short 
distances. On the other hand, the predicted growth of 
the fiux-tube diameter is manifest in the lattice data. 



V. CONCLUSION 

The gluonic-distribution inside the meson has been 
revealed. The Monte- Carlo simulations have been per- 
formed on SU(3) gauge-group for temperatures T ~ 
0.8 Tc and T ~ 0.9 Tc. Noise reduction has been 
achieved by a gauge-independent high statistics ap- 
proach, in addition to the employment of adequate levels 
of gauge smoothing that preserve the relevant physics at 
large distances. This method is variant to noise-reduction 
by Abelian gauge-fixing. The flux tube, characterized as 
a correlation between the action-density and the mesonic 
operator (Polyakov-lines) , has been displayed up to dis- 
tances of 1.4 fm. The flux iso-lines and iso-surfaces 
display a curved profile along the tube. The profile 
is showing a non-uniform action-density pattern unlike 
that observed using Wilson's loop as a mesonic operator 
at T = 0. The flux tube width proflle is compared to 
the corresponding mean-square width of the free Bosonic 
string fluctuations at all planes between the color sources. 
For source separation distances R > 0.7 fm, measure- 



ments of the tube cross-section at each selected trans- 
verse plane show a non-constant width for the tube with 
maximum width at the symmetry point of the tube. At 
small source separation 0.7 fm < R < 1.1 fm, the tube 
is seen to yield a more compact (squeezed) form than the 
string model would predict. The deviations of the tube 
width proflle from the corresponding string proflle de- 
crease gradually as source separation increases and the 
proflles begin to compare well at i? i2 1.3 fm. The glu- 
onic profiles displayed in this investigation are geometric 
manifestations of thermal effects on the qq potential (the 
measured decrease in the string tension). Moreover, the 
squeezed gluonic profile in comparison to the rounded 
string fluctuations provides also a geometrical interpre- 
tation for the deviations of the predicted string tension 
based on the free string picture from the corresponding 
lattice results [l[. 

This study is motivating further investigations of the 
energy-density and chromo-electromagnetic distributions 
with methodological improvements that minimize the 
number of smearing sweeps and increase the number of 
measurements. It would also be interesting to confront 
these profiles with the Bosonic string profiles in the con- 
text of string's self-interactions. The string's geometrical 
effects (curved profiles) ought to be addressed in other 
gauge groups. Work is progressing in these directions, 
in addition to the detailed investigation of the thermal 
hadronic gluonic distributions by straight forward gener- 
alizations of the calculations presented here to Baryons 
;35, 36]. 
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